home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / os2 / adaptor.zip / ADAPT.ZIP / adaptor / examples / laplace / laplace.f < prev    next >
Text File  |  1993-03-23  |  1KB  |  50 lines
  1.       PROGRAM LAPLACE
  2.       REAL F(:,:), DF(:,:)
  3.       REAL HF(:,:)
  4.       LOGICAL CMASK (:,:)
  5.       INTEGER MAXX, MAXY
  6.       REAL FMAX
  7.       INTEGER ITER
  8. c 
  9. c     read in sizes 
  10. c
  11.       PRINT *,'MAXX = (z.B. 64) '
  12.       READ *,MAXX
  13.       PRINT *,'MAXY = (z.B. 64) '
  14.       READ *,MAXY
  15.       ALLOCATE (F(MAXX,MAXY), DF(1:MAXX,MAXY))
  16.       ALLOCATE (HF(MAXX,1:MAXY))
  17.       ALLOCATE (CMASK(MAXX,1:MAXY))
  18. C
  19.       CMASK = .FALSE.
  20.       CMASK (2:MAXX-1,2:MAXY-1) = .TRUE.
  21.       F = 2.
  22.       F(:,MAXY) = 1.
  23.       WHERE (CMASK)
  24.           F = 0.0
  25.       ENDWHERE
  26.       ITER = 0
  27.       FMAX = 1
  28.       DO WHILE (FMAX .gt. 0.001)
  29.          ITER = ITER + 1
  30.          HF = CSHIFT (F,1,1)
  31.          DF = CSHIFT (F,1,-1)
  32.          HF = HF + DF
  33.          DF = CSHIFT (F,2,1)
  34.          HF = HF + DF
  35.          DF = CSHIFT (F,2,-1)
  36.          HF = (HF + DF) * 0.25
  37.          DF = 0
  38.          WHERE (CMASK)
  39.             DF = HF - F
  40.             F  = F + DF
  41.          ENDWHERE
  42.          DF = ABS(DF)
  43.          FMAX = MAXVAL (DF)
  44.          PRINT *,'Iteration ',ITER,'  Max = ',FMAX
  45.       END DO
  46.       PRINT *, ITER, ' Iterationen benoetigt'
  47.       DEALLOCATE (CMASK, HF, DF, F)
  48.       END 
  49.  
  50.